The Hoist Spin Foam Model 
via Cubulations 

A. Baratin 1 * C. Flori 1 J T. Thiemann 1 ' 2 * 

1 MPI f. Gravitationsphysik, Albert-Einstein-I nstitut, 
Am Muhlenberg 1, 14476 Potsdam, Germany 

2 Perimeter Institute for Theoretical Physics, 
31 Caroline Street N, Waterloo, ON N2L 2Y5, Canada 

Preprint AEI-2008-093 



Abstract 

Spin Foam Models (SFM) are an attempt at a covariant or path integral formulation of canonical Loop 
Quantum Gravity (LQG). Traditionally, SFM rely on 1. the Plebanski formulation of GR as a constrained BF 
Theory. 2. simplicial triangulations as a UV regulator and 3. a sum over all triangulations via group field 
techniques (GFT) in order to get rid off triangulation dependence. 

Subtle tasks for current SFM are to establish 1. the correct quantum implementation of Plebanski's 
constraints. 2. the existence of a semiclassical sector implementing additional Regge constraints arising from 
simplicial triangulations and 3. the physical inner product of LQG via GFT. 

We propose a new approach which deals with these issues as follows: 1. The simplicity constraints are 
correctly implemented by starting directly from the Hoist action which is also a proper starting point for 
canonical LQG. 2. Cubulations are chosen rather than triangulations as a regulator. 3. We give a direct 
interpretation of our spin foam model as a generating functional of n - point functions on the physical Hilbert 
space at finite regulator. 

This paper focuses on ideas and tasks to be performed before the model can be taken seriously, however, it 
transpires that 1. this model's amplitudes differ from those of current SFM, 2. tetrad n-point functions reveal 
a Wick like structure and 3. the restriction to simple representations does not occur automatically but must 
employ the time gauge just as in the classical theory. 
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1 Introduction 



Spin foam models (SFM) [I] are an attempt at a covariant or path integral formulation of canonical Loop 
Quantum Gravity (LQG) [21 [Hid]. In their current formulation, SFM exploit the Plebanski formulation [5] 
of pure General Relativity (GR) as a constrained BF theory. This approach is well motivated because one 
can view the Plebanski action as a kind of perturbation of the BF action (albeit the perturbation parameter 
is a Lagrange multiplier field which one needs to integrate over in a path integral). The path integral for 
BF theory, however, is under good control [6j so that one may hope to get a valid path integral formulation 
for GR by functional current derivation methods [7] familiar from ordinary QFT. 

As we will try to explain in the next section (see also [3] ) the quantum implementation of the so called 
simplicity constraints of Plebanski theory, to the best knowledge of the authors, has still not been achieved 
to full satisfaction from first principles in SFM They are called simplicity constraints because they enforce 
the B field of BF theory to be simple, that is, to originate from a tetrad. Clearly, unless the simplicity 
constraints are properly implemented, the resulting theory has little to do with quantum gravity. An issue 
to keep in mind is that the solutions to the classical simplicity constraint consist of five sectors, two of 
which give rise to ± times the Palatini action, two of which give rise to ± times a topological action and a 
degenerate sector. All of these sectors are included in a sum over Plebanski histories which may or may not 
be what one wants. 

It is appropriate to mention also further constraints in SFM at this point: Namely, as we also explain 
in the next section, SFM rely on a simplicial triangulation r of the differential 4-manifold as well as a 
dual graph r*. A recent analysis has shown [8] that freely specifying geometrical data (areas or fluxes) on 
the faces of r tends to lead to inconsistent values of the lengths of the edges of r unless so called Regge 
constraints in addition to the simplicity constraints are imposed. These constraints are important to be 
taken care of if one wants to relate SFM to the established theory of Regge calculus [9] and in order to 
capture the correct semiclassical limit. The underlying reason for these constraints is that Regge calculus is 
formulated directly in terms of edge lengths while in SFM one rather works with electrical fluxes or areas. 
However, a typical simplicial triangulation has far more faces than edges in r so that assigning a length to 
an edge from given area values maybe ambiguous and/or inconsistent. The authors of [8] speculate that this 
may pose a problem also for canonical LQG and that the Hamiltonian constraint of canonical LQG should 
maybe commute withe these Regge constraints. In order avoid further confusion in the literature on this 
puzzling issue we will take here the opportunity to show that there are no Regge like constraints in canonical 
LQG to be taken care of. Notice that since our path integral is explicitly based on the Hoist action, there 
is no necessity to relate it to the Regge action which for Plebanski's theory is of course a challenge. 

One possibility to make progress on these issues is based on a very simple idea that, to the best knowledge 
of the authors, occurs for the first time in [TO]: Namely, simply try to formulate the path integral in terms 
of the Hoist action [11] rather than the Plebanski action. Not only is the Hoist action a valid starting 
point for canonical LQG, also the simplicity constraints are explicitly solved in that one works entirely 
with tetrads from the beginning. More precisely, the Hoist action uses a specific qudratic expression in the 
tetrads for the B field of BF theory which also depends on the Immirzi parameter |25j . Hence, the Hoist 
action depends on a specific, non degenerate linear combination of the four non degenerate solutions of the 
simplicity constraints (see next section for details) and is thus at the same time more general and more 
restricted because the Hoist path integral will not sum over the afore mentioned five sectors of Plebanski's 
theory. As already mentioned, it is at present debated how the fact that one actually takes a sum over all 
histories with a mixture of positive and negative Palatini and topological actions affects the semiclassical 
properties of the Plebanski path integral. 

As observed in [TO], since the Hoist action is quadratic in the tetrads, one can in principle integrate 
out the tetrad in the resulting Gaussian integral. This has been sketched in [10] , however, the expressions 
given there are far from rigorous. Here we will give a rigorous expression. Also, we will include the correct 
measure factor |12] resulting from the second class constraints involved in the Hoist action which makes sure 
that the path integral qualifies as a reduced phase space quantisation of the theory as has been stressed 
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in [13] . A similar analysis has been carried out for the Plebanski theory in [13], however, the resulting 
measure factor is widely ignored in the SFM literature. The result of the Gaussian integral is an interesting 
determinant that displays the full non linearity of Einstein's theory. When translating the remaining integral 
over the connection in the partition function into SFM language, that is, sums over vertex, edge and face 
representations, one sees that our model differs drastically from all current SFM. 

Of course, we also need to introduce an IR and UV regulator in the form of a finite cell decomposition. 
Two observations lead us to depart from the usual SFM approach where one works with simplicial cell 
complexes. The first is the result [15] which demonstrates that current semiclassical states used in LQG 
do not assign good classical behaviour to the volume operator [16] of LQG unless the underlying graph has 
cubic topology. Since the volume operator plays a pivotal role for LQG as it defines triad operators and 
hence the dynamics, this is a first motivation to consider cubic triangulations of the four manifold, also 
called "cubulations" . Notice that any four manifold can be cubulated and that within each chart of an atlas 
the cubulation can be chosen to be regular (see e.g. [16] and references therein). The second observation 
is that the original motivation for considering simplicial cell complexes in current SFM comes from their 
closeness to BF theory. BF theory is a topological QFT and therefore one would like to keep triangulation 
independence of the BF SFM amplitude. That this is actually true is a celebrated result in BF theory. 
In particular, in order to keep triangulation independence it is necessary to integrate the B field over the 
triangles t of the tringulation and the F field over the faces / bounding the loops in a dual graph [19] . 
However, GR is not a TQFT and therefore the requirement to have triangulation independence is somewhat 
obscure. Of course it is natural if one wants to exploit the properties of BF theory but not if one takes a 
different route as we tend to do here. Hence, if we drop that requirement, then it is much more natural to 
refrain from considering the dual graph in addition to the triangulation. This also solves another issue: 
Notice that the gauge group SO(p,q) acts on the B field of BF theory by the adjoint action and on the 
connection A underlying F in the usual way. The question is where the gauge transformation acts on the 
discretised variables B(t), A(df) (flux and holonomy). It would be natural to have the gauge group act 
at the barycentres of t and at the starting point of the loop df which will be a vertex of the dual graph. 
However, notice that the vertices of the dual graph and the triangles are disjoint from each other, the edges 
of the graph are dual to the tetrahedra of the cell complex. Hence, local gauge invariance in discretised BF 
theory at the level of the action is not manifest and even less in Plebanski theory. In fact gauge invariance 
is related to the closure constraint in SFM which is a subtle issue as we will see in the next section. If one 
works just with a triangulation and drops the dual graph then gauge invariance issues are easy to take care 
of. Hence, it is motivated to work with a triangulation that maximally simplifies the Gaussian integral. As 
we will show, this again leads to cubulations. This also nicely fits with the framework of Algebraic Quantum 
Gravity [20] which in its minimal version also is formulated in terms of algebraic graphs of cubic topology 
only. 

The architecture of this article is as follows: 

In section two we give a non technical review of current spin foam models. We sketch their derivation 
from the classical Plebanski action focussing on the points where a first principle argument is missing. These 
issues will be the motivation for our different route. 

In section three we clarify the issue mentioned above, namely whether the Regge constraints on SFM 
discovered recently in [8] also arise in the canonical formulation. We prove that they do not play any role 
in LQG. 

In section four we derive the Hoist SFM using cubulations as UV regulator as motivated above. As this 
is a exploratory paper only, we will not worry about convergence issues which will be properly addressed 
in subsequent works. More precisely, what we compute are tetrad n-point functions. These should contain 
sufficient information to compute anything of interest in LQG such as graviton scattering amplitudes as in 
|21j via LSZ (Lehmann - Zimmermann - Symanzik) like formulas as in ordinary QFT [22j which allows 
to reconstruct the S - matrix from symmetric vacuum n - point functions. Of course, how these n-point 
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functions are related to true observables in a diffeomorphism invariant theory is a subtle issue which will be 
clarified in a separate paper [23] ■ Here we only give a summary. The n-point functions can be computed 
in closed form up to a remaining functional integral over the connections. This can be done for either 
signature of the spacetime metric. At this point one could invoke SFM techniques and expand the integral 
using harmonic analysis on the gauge group. The resulting intertwiner displays a much more complictated 
structure than in any of the current spin foam models. In particular, pictorially speaking one basic building 
block is an octagon diagramme an analytic expression for which could be called the 96 - j symbol in the 
case of G = SO{A). Yet, the n-point functions display a certain Wick like structure as if they came from a 
Gaussian integral. What makes the theory interacting and obstructs the tetrads from being a generalised 
free fielcQ is the additional functional integral over the connection. In background dependent QFT the 
moments of a Gaussian measure depend on a background dependent covariance (usually depending on the 
Laplacian (in the Euclidean setting) and the mass). Our theory behaves similar, just that due to background 
independence the covariance is itsef a field that must be integrated over. This is similar in spirit to what 
happens in 3D [24] when coupling GR to point particles: There, when integrating over the gravitational 
degrees of freedom one ends up with particles moving on a non commutative geometry. Here instead of a 
non commutative geometry we obtain an interacting theory of tetrad fields. 

In section five we conclude and outline the missing tasks that need to be performed before our model can 
be taken seriously. An interesting result of our analysis is that in the present formulation which lacks the 
simplicity constraints of Plebansk's theory, the irreducible representations of Spin(p, 4 — p) are not forced 
to be simple. Simple representations, which basically reduce Spin(p,4 — p) to an SU(2) subgroup, can only 
arise if we impose the time gauge which in the classical theory is used in order to reduce the Hoist connection 
to the Ashtekar - Barbero - Immirz connection which is a necessary ingredient in the canonical quantisation 
programme. Gauge fixing conditions of course naturally arise in any attempt to make formal path integral 
expressions better behaved and here the situation is similar. 

Two appendices treat some simple technical aspects of this work. 

Most parts of this paper do not depend on whether the spacetime signature is Lorentzian or Euclidean. 

2 Outline of Current Spin Foam Models 

In this section we intend to give a brief summary of the developments in SFM theory with a focus on the 
derivation of the current models from the Plebanski action and the gaps in that derivation. This serves as 
the motivation for the present paper. 

To begin with, it is worth mentioning that the classical solutions to the simplicity constraints actually 
comprise altogether five sectors, namely two topological sectors B = ±e A e, two Palatini sectors B = *e A e 
(where * denotes the Hodge map with respect to the internal Minkowski or Euclidean metric) and one 
degenerate sector. In a path integral a sum over all these sectors will occur while one would expect that one 
should only include one of the Palatini sectors or maybe a Hoist sector B = *eAe + ^-eAe [11] in order to 
have a path integral for Einstein's theory. Here 7 is the Immirzi parameter of LQG [25] . 




where $ is a scalar Lagrange multiplier field with values in the tensor product of two copies of so(l,3) or 
so(4) depending on the signature and F is the curvature of the connection A. In a formal path integral 
formulation one integrates exp(iS) over A, B, $. Integrating first over we are left with a partition function 
of the form 




We now sketch the usual "derivation" of SFM from the Plebanski action: 
The Plebanski action is of the form 







1 Roughly, a generalised free scalar field is such that all its n-point functions are already determined by its two point function. 
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where C(B) denotes the collection of the simplicity constraints on B. If one would solve the delta distribution 
by integration over B one would get the afore mentioned sum over the five sectors and integral over the 
tetrad fields. However, this would result in a complicated expression which does not exploit the relation of 
Plebanski's formulation to BF theory. Thus, rather than doing that, one notices that roughly speaking 

Bexp(ij Tr(B A F)) = i — exp(i J Tr(B A F)) (2.3) 

Denoting the functional derivative by X one can now formally pull the 5 distribution out of the B integral 
and perform the integration over B resulting in 

Z = j [dA] 5(C(X)) ■ 5(F) (2.4) 

Without the "operator" 5(C(X)) this would be the formal partition function of BF theory. Thus, one has 
achieved the goal to preserve the closeness of the theory to BF theory. One now should expand 5(F) in 
terms of eigenfunctions of the collection of operators C(X) and keep only the zero eigenfunctions multiplied 
by<5(0). 

In order to give meaning to those formal expressions one has to introduce a UV and IR cutoff as is 
customary in constructive QFT. That is, one considers finite simplicial triangulations r of the (possibly 
compact) differential 4 manifold and dual graphs r*. The two forms B are now approximated by integrals 
B(t) of B over triangles t of r while the curvatures F are approximated by holonomies A(df) around the 
loops df of the faces / dual to the triangles t. One writes f(t) for the face dual to t. The BF action is then 
discretised by 

£ Tr(A(df(t)) B(t)) (2.5) 
The reason to work with both r and r* is that in fact 

/ Tr(B AF) = V Tr(F(f(t))B(t)) (2.6) 
Jm ter 

is an exact identity [19] where F(f) denotes the integral of F over /. This is very convenient in particular 
for pure BF theory. The only approximation thus consists in replacing F(f) by A(df) — \q- 

Likewise, the functional derivatives X must be approximated by ordinary derivatives Xt with respect 
to the variables F(t) := A(df(t)) — A(df(t)) T . Notice that when defined like that, the Xt are mutually 
commuting. The exploration of the model with this definition of Xt will appear in [26]. However, this is not 
what is done in current SFM. Rather one replaces Xt by Yt, the right invariant vector field on the copy of 
G associated with the variable A(df(t)). The reason for doing that is that the Yt have a simpler action on 
the delta distribution 

5(F):=Yl5 G (A(df(t)) (2.7) 

It is usually justified by saying that 5(F) has support on A(df) = 1q and that Y t , X t differ by multiplication 
with holonomies which should be supported at 1q- However, this argument is certainly not rigorous because 
the support of the 5 distribution can drastically change when acting with differential operators. Moreover, 
as already said, this substitution comes with a price: While the simplicity constraints in terms of Xt 
are mutually commuting, those in terms of Yt are not, in fact they are inconsistent with each other. In 
addition, one does not impose all the simplicity constraints but only a subset of them: There are three 
types: Constraints involving 1. the same triangle, 2. two triangles sharing an edge and 3. two triangles 
sharing a vertex. The latter constraint is implied by the so called closure constraint on tetrahedra T 

£ Y t = (2.8) 

ter 
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(but not vice versa). This constraint looks as if it would be automatically satisfied because it looks like a 
gauge invariance condition. However, the product of 5 distributions (|2,7p in A(i ? ) is not annihilated by the 
closure constraint (|2.8p ! This is obvious from the fact that the product of 5 distributions involves products 
of the form 

Hx,M(df(t)) (2.9) 

teT 

where ir denotes an irreducuible representation of G and Xn its character. However, there is no gauge 
invariant interwiner among the loops df(t), t E T. One usually argues that the closure constraint is taken 
into account because after integrating over A one is only left with gauge invariant intertwiners, but this is 
wrong before integrating. In fact since integration with the Haar measure always projects out the gauge 
invariant part, anything can be made gauge invariant this way. We feel that neglecting the 3rd kind of 
simplicity constrant (implied by taking the closure constraint for granted) makes the model too local. The 
effect of truly taking the closure constraint into account will also be explored in |26j. 

As already said, even the simplicity constraints of the first two types are anomalous as they imply 
vanishing volume [3] and fix the above mentioned intertwiner to be unique (the model has not enough 
degrees of freedom). This and other investigations has ruled out the Barrett - Crane model [6] which 
however was an important step in SFM research because it triggered the development of model independent 
mathematical tools. Recent activities in SFM therefore have focussed on trying to implement the simplicity 
constraints of the first two types differently. There are two camps: In [27] one uses Master Constraint type of 
techniques which were developed in a different context [28J in order to treat anomalous (rather: second class) 
constraints. In [29] one refrains from imposing the simplicity constraints as operator conditions altogether 
but rather imposes them semiclassically by expanding SFM amplitudes in terms of group coherent states 
|30j developed by Perelomov [31J and then uses the simplicity condition on the classical bivectors on which 
the semiclassical amplitudes depend. 

In our view, while these methods seem to give rise to models with better semiclassical properties |32] and 
to a better understanding whether one is looking at the topological or Palatini sector respectively, in our 
opinion it is fair to say that a satisfactory derivation from first principles is missing. By this we mean that 
one should be able to arrive at those models starting from the Plebanski action, another classically equivalent 
action or the Hamiltonian formulation and then carry out integrations and imposition of constraints without 
intermediate approximations or ad hoc substitutions as those listed above. This is precisly the motivation 
of the current paper. 



3 On Regge Constraints in Canonical LQG 

In canonical LQG the states are not labelled by a pair consisting of a graph and a dual polyhedronal complex 
(in fact, generic graphs do not admit such a dual complex at all). Rather states are just labelled by graphs, 
edge representations and vertex intertwiners. This is the spin network orthonormal basis. In fact, all graphs 
are used in canonical LQG, not only four valent ones as the ones that are used in present SFM analysis based 
on simplicial trangulations which is why current SFM only capture a small subspace of the LQG Hilbert 
space. This subspace is not left invariant under the LQG holonomy flux algebra [33j so that this subspace 
is not a representation of it. Hence, in order to establish the relation between LQG and SFM it would seem 
that one needs to include spin network states on all possible boundary graphs into the SFM analysi^H. 

In addition to this basis of states, one has holonomy, flux, area, triad or length operators labelled by 
curves or surfaces. None of them has a physical meaning, it is only their occurence in compound operators 
assembeled from them and which are Dirac observables or constraint operators respectively which has a 
physical meaning. These arise as integrals over spatial scalar densities and by the theorem established in 
paper V of the series [35] any dependence of these quantities on the electric degrees of freedom is via triad 
like operators. More in detail, the cylindrical projections of such operators are labelled by arbitrary graphs 

2 Except if one follows the philosophy of |20| . 
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and the triad like operators as well as the holonomies that arise in the corresponding expressions are defined 
by the edges of that graph. Hence we see that also physically relevant operators just depend on the structure 
of a graph and not on a dual cell complex. 

This is maybe confusing because flux and operators depending on arbitrary (piecewise analytic) surfaces 
and curves actually define the fundamental kinematical algebra underlying LQG. However, flux operators 
are used at intermediate stages, they are not even Gauss invariant. What one uses are volume operators 
and constructs from them length operators [M] and more generally triad like operators via the techniques 
developed in |35j. These triad operators which enter gauge invariant and diffeomorphism co(in)variant 
expressions such as the Hamiltonian constraint only excite the states along the edges of its underlying 
graphs. Again, there is no reference to any polyhedronal cell complex whatsoever. There are efforts to 
invent a length operator which also uses a cell complex [36], however, the resulting operator in contrast to 
|34| is not cylindrically consistent. 

Now although they do not have a direct physical meaning, LQG admits these holonomy and flux or 
area operators labelled by arbitrary curves and surfaces respectively not related to any graph [161 139] . A 
given collection of such operators may or may not have good semiclassical properties with respect to a 
given semiclassical state. But this is always true in quantum theory: Even the standard coherent states of 
the harmonic oscillator do not assign good semiclassical expectation values to all operators that one can 
construct in the corresponding Hilbert space, see e.g. [40] for a detailed discussion. What is important 
is that there is a sufficient number of functions on phase space that should seperate the points^]. In our 
case these could be holonomies and triads integrated along the edges of the graph underlying the given 
semiclassical state. These variables are enough to separate the points of the classical phase space to an 
arbitrary precision depending on how fine the given graph is embedded with respect to the background 
metric to be approximated. 

Now it is true that there are certain semiclassical states in LQG that depends on a polyhedronal 3D 
complex r |37] as an additional structure. However, this 3D r does not have any dynamical interpretation 
as in SFM since it does not define a 4D extension. It just serves to construct a complexifier [37J and thus 
to produce an infinite family of states labelled by r, a fixed point (Aq, Eq) in the classical phase space and 
any graph 7. Any sufficiently fine r would do and so r does not have any fundamental meaning. In fact one 
could choose to consider complexifiers that do not come from any r, see e.g. the first reference of [37] . The 
semiclassical properties are optimal for those members of the family for which 7 is in fact dual to r. For 
each edge e of such a graph there is a classical flux value Eo(S e ) and a classical holonomy Ao(e) which label 
the state where S e is the unique face in r dual to e. Notice that what one is interested in are the lengths 
of the edges e in 7 and not the length of the links in the 1-skeleton of r. Due to duality, there is a one to 
one correspondence between the classical values Eq(S £ ) and the expectation values of the length operators 
associated to e, no overcounting of degrees of freedom takes place, no ambiguities arise, whatever r is. 

That the corresponding length operators assume the correct semiclassical expectation value with respect 
to these states follows from the analysis performed in [381 [20] where one has chosen cubical r for simplicity. In 
fact, the analysis of [15] has demonstrated that the volume operator of LQG assumes the correct semiclassical 
expectation values only if r has cubical topology so that one does not really have too much of a choice as far 
as r is concerned. This in turn implies that current spin foam models based on simplicial cell complexes do 
not admit the semiclassical states [37j as boundary states which could mean that the current models maybe 
have to be extended to more general triangulations. 

Notice that curiously, in 3d, in contrast to simplicial cell complexes, the number of faces of (sufficiently 
regular) cubical r equals the number of links in the 1-skeleton of r so that also here there is a one to one 
correspondence between faces and links within r. However, this does not have any relevance for our main 
point, namely that in LQG the fundamental quantum objects are holonomies along the edges of graphs and 

3 In LQG all semiclassical states depend on a (possibly infinite) graph due to the non separability of the Hilbert space. 
Therefore one cannot require good semiclasical bahaviour for all holonomies but only for those associated to the given graph. 
That graph should be sufficiently densely embedded in order that one can reconstruct the connection from the given holonomies 
up to an error that depends on the fraction of the graph scale divided by the curvature scale, see [37] for an extensive discussion. 
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triad or length operators associated to those edges, there is no association of length or area operators to 
links or faces in some dual cell complex. 
To summarise: 

In LQG there is a basis of states labelled by graphs. Given an arbitrarily small error, there is a set of operators 
whose classical counterparts separate the points of the classical phase space up to that error and whose 
cylindrical projections are labelled by the edges of graphs. In addition, there are also semiclassical states 
depending on graphs with respect to which this separating system of operators have optimaQ semiclassical 
properties. There is no triangulation with respect to which one would be interested in relating the lengths 
of the edges of its 1 - skeleton to the areas of its surfaces. Regge like constraints never occur. 

4 Derivation of the Model 

This is the main section of the present paper. It is subdivided into five parts. In the first we motivate the 
use of cubulations from different perspectives and discuss some of their properties. In the second we sketch 
the relation between path integral n-point functions and physical (observable) correlators in terms of the 
physical inner product of the theory. More details on that issue will appear elsewhere |23| . This crucially 
works via a choice of gauge fixing or clock system. In the third part we apply our machinery to n-point 
tetrad functions or equivalently to a generating function of a (complex, regulated) measure. That measure 
displays a Gaussian like structure and we can accordingly integrate out half of the degrees of freedom under 
some assumptions about the choice of gauge fixing. In the fourth part we discuss the properties of the 
resulting integral over the remaining degrees of freedom, its Wick like structure and the structure of the 
vertex amplitude of the corresponding SFM model when invoking harmonic analysis on the gauge group. 
Finally, in the fifth part we discuss how these n-point functions are related to the physical inner product and 
the kinematical Hilbert space of LQG, in particular, how the covariant connection of the Hoist path integral 
reduces to the Ashtekar - Barbero - Immirzi connection of the canonical theory in physical amplitudes. 

4.1 Cubulations 

We start by motivating the use of cubulations rather than simplicial triangulations. 
4.1.1 Gauge invariance 

Let us look more closely at the issue of gauge invariance for BF theory which makes use also of a dual 
graph. Here gauge invariance is not preserved locally (i.e. triangle wise) in the formula J Tr(B A F) = 
^ t Tr(i?(t)i ? (/(t)) if both B and F transform locally in the adjoint representation. In order to make the 
gauge transformations more local, one could discretise them. To see how this can be achieved, recall that by 
definition of a cell dual to a simplesEI in a simplicial complex r, the face f(t) is a union of triangles [i, T, a] 
subject to the condition t C dT, T C da. Here (.) denotes the barycentre [19] of a simplex and T, a denote 
tetrahedra and four simplicies in r respectively. So we see that both t and f(t) contain the barycentre t 
in their intersection and we could define a disjoint action of the gauge group on both B(t), F(f(t)) at t. 
However, this is no longer possible when using the approximation ^ t Tr(i?(t)A(3/(t))) because now the 
only natural action of the gauge group on the loop holonomy is by adjoint action at a starting point on 
df(t). Now df(t) is a composition of the half edges [T,a] where t C dT, T C da but the fundamental 
degrees of freedom are the holonomies along the edges e = [a, a'] for a Da' = T, t C dT. Obviously, the only 
natural starting point of the loops is then at the vertices a which are disjoint from the triangles t. But the 
triangles are also disjoint from the half edges as a simple calculation reveals. Hence, in order to maintain 
gauge invariance one has to invent an unnatural discretised action of the gauge group. 

4 In the sense of saturating the Heisenberg uncertainty relation. 
As usual |19] , an n-simplex is denoted by [po, --,Pn] where the points Pi denote its corners. 
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We do not know if such a consistent prescription can be found at all but rather take these complications 
that come from the fact that one is dealing simultaneously with a (simplicial) complex and its dual cell 
complex as a further piece of motivation to work just with the triangulation. 

4.1.2 Cubulations versus simplicial triangulations 

The previous considerations do not specify the type of triangulations to be considered. As already said, 
the first piece of information why to use cubulations rather than simplicial triangulations is because the 
boundary graphs must contain cubical ones in order to make sure that the corresponding boundary Hilbert 
space contains enough semiclassical states [T5]. However, there is an additional, more practical motivation 
for doing so which we discuss now. 

Recall that the Hoist action is given by 

S = — f Tr(G[A] A e A e) = - / G U [A] A e 1 A e J (4.1) 

Here k denotes Newton's constant, 

G[A] = 2(*F[A] + -F[A\) (4.2) 

L " 7 

where Fjj = dAj j + Ajk A A K j denotes the curvature of the connection A, 7 is the Immirzi parameter, * 
denotes the internal Hodge dual, that is, 

(*T)jj := ^e IJKL r, KM V LN T MN (4.3) 

where I,J,K,.. = 0, ..,3 and 77 is the Minkowski or Euclidian metric for structure group G = 50(1,3) or 
G = SO (4) respectively. As motivated in the introduction, we plan to keep the co - tetrad one forms e 1 
rather than introducing a B field and thus the simplicity constraints are manifestly solved. Moreover, the 
issue raised in [8] is circumvented as co-tetrads are labelled by curves and not by (overcomplete) surfaces. 

In order to give meaning to a path integral formulation we consider a UV cutoff in terms of a triangulation 
r of M which we choose to be finite, thereby introducing an IR regulator as well. Let us denote the two 
- dimensional faces of r by / and the one dimensional edges of r by I. We want to discretise (|4.ip in a 
manifestly (and locally) gauge invariant way, just using edges and faces. To do so we equip all edges with 
an orientation once and for all. Given an edge I consider 

4 := J [A(!(x))? je J (x) (4.4) 

Here l(x) for x £ I denotes the segment of I that starts at the starting point of I and ends at x and 
L4(p)] J denotes the G valued holonomy of A along a path p. Evidently, under local gauge transformations 
g : M —> G, (|4.4p transforms as e\ \—* g j(b(l)) ef where b(l) denotes the beginning point of I. 

To avoid confusion, here g G G means the following: The fundamental objects are the matrices g 1 j. 
Set gu := r] IK g K j. Then g G G iff giK~gjm KL = mJ- This is equivalent with (g^ 1 ) 1 j = r] IL g K l i]kj- In 
other words 

(g- 1 ) = (~gf (4-5) 

If g 1 j = [e~x.p(F)Y j for some generator F 1 j then (|4.5p means that Fjj + Fjj = 0. In abuse of notation 
one usually uses the same symbols g, F and g, F respectively but unless we are in the Euclidian regime we 
should pay attention to the index position. 
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Clearly, the curvature F must be discretised in terms of the holonomy of A along the closed loops df 
where we have also equipped the faces / with an orientation once and for all. We have 

Fu(f) := ^A^l/J-PiS))- 1 ]^) 

= \{[Mdf)\u - [{(Mff)) T ]u) 




if 



(4.6) 

where we have used the non Abelian Stokes theorem for "small" loops, that is 

A(df) « exp( / F) (4.7) 



Jf 

and we have written Fu{x) := Fjj(x). We may now define the antisymmetric matrix 

Gij(f) = (*F(f))u + -Fu(f) (4.8) 

7 

Imagine now that we would use a simplicial triangulation. Hence M is a disjoint (up to common 
tetrahedra) union of four simplicices o = \po(a), ..,^4(0")]. For each pj(a) label the four boundary edges of a 
starting at Pj(cr) by lji(cr) and let the face (triangle) of a spanned by ljj,(cr), ll(cr) be denoted by fjiv(<j) with 
the convention /^(cr) = —fl^ia)- Now the orientation of Zju(cr) either coincides with the given orientation 
of the corresponding edge in a or it does not. In the former case define ejf(a) := eL while in the latter 

we define ejf(a) := [A(lf l (a))~ 1 e lJ '^Y ■ Then we have 



kS = [ Tr ( G A e A e ) 

<t€t 

« \ E E ^ pX Gu{fl u {a))e^{a) e{\a) 



=: J]G/^ (4.9) 

i,l' 

where we have averaged over the corners of a 4 - simplex. For any simplicial triangulation the (symmetric 
in the compound index (1,1)) matrix G l jj is difficult to present explicitly due to bookkeeping problems, 
even if we refrain from averaging over the five corners of a 4 - simplex. Moreover, as we intend to perform 
a Gaussian integral over the e\ , we need the determinant of that matrix which is impossible to compute 
explicitly unless it is block diagonal in some sense. 

The latter observation points to a possible solution. First of all any manifold admits a cubulation, that 
is, a triangulation by embedded hypercubea^ [17]. We assume that M has a countable cover by open sets 
O a . Consider a stratification by 4D regions S a subordinate to it. Then S a admits a regular cubulation, 
that is, the 1 - skeleton of the cubulation of M restricted to S a can be chosen to be a regular cubic lattice. 
Non trivial departures from the regular cubulation only appear at the boundaries of the S a . We restrict 
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An easy proof uses the fact that every manifold can be triangulated by simplices. Given a D - simplex, consider the 
barycentre of each of its (^±) sub - p - simplices for p — 0, .., D. Connect the barycentre of any p + 1 - simplex with the 
barycentres of the p - simplices in its boundary. It is not difficult to see that this defines a cubulation of the D - simplex and 
that all p - cubes thus defined are the same ones in common q - simplices of the original simplicial complex. In other words, 
every simplicial complex has a cubulated refinement. 
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attention to those M admitting a cubulation such that in every compact submanifold the ratio of the number 
of cubes involved in the non - regular regions divided by the number of cubes involved in the regular regions 
converges to zero as we take the cubulation to the continuum. For those M, up to corrections which vanish 
in the continuum limit we can treat M as if it would admit a global, regular cubulation. 

Given a regular cubulation r, consider its set of vertices. In 4D, each vertex v is eight valent and there 
are four pairs of edges such that the members of each pair are analytic continuations of each other while 
the tangents at v of four members from mutualy different pairs are linearly independent of each other. It 
is therefore possible to assign to each edge a direction fj, = 0,1,2,3 and an orientation such that adjacent 
edges in the same direction have a common analytic continuation and agree in their orientation. We label 
the edges starting at v in fj, direction by l^(v). Notice that this labelling exhausts all possible edges and 
unambiguously assigns an orientation to all of them. The discretised co - tetrad is then given by 

e» == (4-10) 

Notice that the hypercubic lattice that results solves all our bookkeeping problems since we now may label 
each vertex by a point in Z 4 . 

Next, given a vertex v we denote by v ± ft the next neighbour vertex in fi direction. We define the 
plaquette loop in the fi, v plane at v by 

dU„(v) := Z» o l v (v + ft) o l^v + z))- 1 o l v ( v y l (4.11) 

so that df„^(v) = [<9/^^(w)] _1 . Notice that again this labelling exhausts all minimal loops in the one skeleton 
of t. The discretised "curvature" is therefore 

GTj(v) := e^G u (f pa (v)) (4.12) 

Denoting by a the 4D hypercubes in r we notice that there is a one to one correspondence between the 
vertices v in the - skeleton of r and the hypercubes given by assigning to a that corner v = (zq, .., z^) of 
a with smallest values of all zq, .., z^ € Z. We then find 



kS = V / GuAe 1 A 
„ Jm 



e J 



-EE G TM 4(v) (4.13) 

v I,J,fj,,u 

The crucial observation is now the following: Assemble pairs of indices into a joint index A = (J, //), B = 
(J, v) etc. and let e (v) := e^(u), Gab(v) '■= Gjj(v) etc.. Notice that by construction Gab(v) = Gba(v) 
for all v. Then (|4.13p can be written as 

This means that using (regular) cubulations indeed the matrix G l \j becomes block diagonal where each 
block is labelled by a vertex and corresponds to the symmetric 16 x 16 matrix G(v). This is what makes 
the computation of the detrminant of the huge matrix with entries G l \j practically possible. As we will see, 
the matrices G(v) have a lot of intriguing symmetries which makes the computation of their determinant 
an interesting task. 

Interesting questions that arise in algebraic topology and which we intend to address in future publica- 
tions are: 

1. Given any D - cubulation, does there exist a cubulated refinement such that one can consistently 
assign to every D cube a a vertex v and to all edges an orientation such that there are precisely D 
edges outgoing from v? We call cubulations for which this is possible regular. If that would be the 
case, we could generalise our discretisation from regular hypercubic lattices to arbitrary cubic ones and 
thus would not have to make any error at the boundaries of the stratified regions mentioned above. 
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2. If the answer to [1.] is negative, can one choose maximally regular cubulations as to minimize the 
error in our assumption of globally regular cubulations? In 3D some results on that issue seem to exist 

3. Given maximally regular cubulations, can one make an error estimate resulting from the neglection of 
the non - trivial topology? 

4.2 Notes on n-point functions 

In the SFM literature the first task that one addresses is the computation of the partition function. However, 
the partition function itself has no obvious physical meaning even if one imposes boundary conditions on 
the paths (spin foams) to be integrated (summed) over. The hope is that SFM provide a formula for the 
physical inner product of the underlying constrained canonical theory which starts from some kinematical 
Hilbert space Tt. The purpose of this section is to sketch the connection between path integrals and n - point 
functions for a general constrained theory. We will use reduced phase space quantisation as our starting 
point. The connection with operator constraint quantisation and group averaging [42J and more details will 
appear elsewhere [23] . 

We assume that we are given a classical theory with first class constraints {F} and possibly second class 
constraints {S}. We turn the system into a purely second class system by supplementing {F} with suitable 
gauge fixing conditions {G}. The canonical Hamiltonian H c is a linear combination of the primary con- 
straints plus a piece H' non - vanishing on the constraint surface of the primary constraints (it could be 
identically zero). It can also be written as a first class piece Hq and (some of) the first class constraints 
F. The gauge fixing conditions fix the Lagrange multipliers involved in the canonical Hamiltonian. One 
may split the complete set of canonical pairs (q,p) on the full phase space into two sets ((f), n), (Q, P) such 
that one can solve the system S = F = G = which defines the constraint surface for (</>, tt) = f(Q, P) in 
terms of Q, P. The Q, P are coordinates on the reduced phase space which is eqipped with the pull - back 
symplectic structural! induced by the embedding of the constraint surface specified by /. 

The gauge fixing conditions also induce a reduced Hamiltonian H r which only depends on Q,P and 
which arises by computing the equations of motion for Q, P with respect to H c and then restricting them to 
the gauge fixed values of the Lagrange multipiers and to the constraint surface. Then H r is defined as the 
function of Q, P onljj§ which generates these same equations of motion. We are now in the situation of an 
ordinary Hamiltonian system equipped with a true Hamiltonian H r . We quantise a suitable subalgebra of 
the reduced Poisson algebra as a *— algebra 21 and represent it on a Hilbert space TL. This Hilbert space is to 
be identified with the physical Hilbert space arising from reduced phase space quantisation. Its relation with 
Dirac's constraint quantisation will appear elsewhere [23J. Let t >— > U(t) be the unitary evolution induced 
by H r . Then the object of interest is the transition amplitude or n-point function 

< ip f , U(t f - t n )a n U(t n - i n _i)a n _i..J7(i 2 - - h)^i > (4.15) 

between initial and final states ipi,ipf at initial and final times tj, tt respectively with intermediate measure- 
ments of the operators a%, .., a n 6 21 at t\ <t2 < ■ ■ < t n . 

Preferrably one would like to be in a situation in which there is a cyclic vector £1 for 21 which is also 
a ground state for H r . The existence of a cyclic vector is no restriction because representations of 21 are 
always direct sums of cyclic representations. In this case 2lf2 is dense in TL and we may therefore restrict 
attention to ipi = tpf = by choosing appropriate a\ , .., a n in (|4,15p . The existence of a vacuum state for 
H r means that zero is in the point spectrum of H r . Let us make this assumption for simplicity. 

Let us abbreviate the Heisenberg time evolution as etfc(i) := U (t)^ 1 akU (t) . In principle it would be 
sufficient to restrict the a& to be configuration operators Q because their time evolution contains sufficient 

7 This symplectic structure coincides with the pull - back of the degenerate symplectic structure on the full phase space 
corresponding to the Dirac bracket induced by the system {S, F, G} [41] . 

8 For simplicity, we are assuming a gauge fixing which leads to a conservative reduced Hamiltonian. 
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information about P as well. However, we will stick to the more general case for reasons that will become 
clear later. This leads us the n-point function 

s &> ">*»>■= <n,u(t f -u)n> (416) 

where we have properly normalised as to give the - point function the value unity. This has the advantage 
that certain infinities that woudl otherwise arise in the following can be absorbed. Notice that since $7 is a 
ground state, the U(tf) and U(U) as well as the denominator could be dropped in (|4.16j) . 

Now a combination of well known heuristic arguments [H] , [33] reviewed in [23] reveals the following: 
Consider any initial and final configuration qi,qf on the full phase space and denote by V((ti,qi), (£/,<?/)) 
the set of paths^l in full configuration space between (%,<?/) at times U,tf respectively. Consider 

Z[j; qi ,q f ]= lim / [Dq Dp DX Dp] S[G] \ det[{F, G}]\ p e i s[q ' pX ^ e* ft t dt m ' q{t) (4.17) 

-**> t /-*°°./p((t i ,, i ),(t/,?/)) 

Here j is a current in the fibre bundle dual to that of q, S[q,p, A, p] is the canonical action after performing 
the singular Legendre transform from the Lagrangian to the Hamiltonian formulation^! and p is a local 
function of q,p which is usually related to the Dirac bracket determinant det[{5, S}] [43J. 

Now the primary constraints are always of the form tt = f(Q,P,(p) where we have split again the 
canonical pairs into two groups. Thus, S[q,p, A, p] is linear in those momenta tt and we can integrate them 
out yielding 5 distributions of the form 6[X — (.)] 5[p — (.)] which can be solved by integrating over X,p. If 
we assume that the dependence of the remaining action on P is only quadratic and that G and | detfji 7 , G}]\ 
are independent of P then we can integrate also over P which yields in general a Jacobian / coming from 
the Legendre transform. We can then write (|4.17p as 

Z[j;qi,q f ]= lim / [Dq] S[G] \ det[{F, G}]\ p I e iS H dt m ' q(t) (4.18) 

-*i.*/-°°77'((t ( ,ci) I (t / ,c / )) 

where proper substitutions of tt from solving the primary constraints and of P from the Legendre transfor- 
mation are understood. Here S[q] is the original (covariant) Lagrangian action. 
Defining x[j] '■= z[o- gff^j ^ e covarian t or path integral n - point functions 

have the canonical or physical interpretation of 

<n,T(a 1 (t l )..a n (t n ))n> (4.20) 

where T is the time ordering symbol, f2 is the aforementioned cyclic vacuum vector defined by the physical 
(or reduced) Hamiltonian H r induced by the gauge fixing G, ak(t) is the Heisenberg operator at time t 
(evolved with respect to H r ) corresponding to and a>k itself classically corresponds to a component of 
q evaluated on the constraint surface S = F = G = 0. The scalar product corresponds to a quantisation 
on the reduced phase space defined by G. Notice how the gauge fixing condition G (or choice of clocks) 
prominently finds its way both into the canonical theory and into the path integral formula (|4.18p . In 
particular, notice that the seemingly similar expression 

Z'[j; qi ,q f ]= lim / [Dq] e if U dt ^>^ (4.21) 

-U,tf->cc J V ((U,q t ),(t f ,q f )) 



9 This should be a suitable measurable space but we leave it unspecified 

10 The Lagrange multipliers A, fi of the primary first and second class constraints respectively play the role of the velocities 
that one could not solve for in terms of the momenta in the process of the Legendre transform 
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does not have any obvious physical interpretation and in addition lacks the important measure factors p, I. 
Remarks: 

1. One may be puzzled by the following: From ordinary gauge theories on background spacetimes such as 
Yang - Mills theory on Minkowski space the path integral or more precisely the generating functional 
of the Schwinger functions (in the Euclidian formulation) does not require any gauge fixing in order 
to give the path integral a physical interpretation. One needs it only in order to divide out the gauge 
volume in a systematic way (Fadeev - Popov identity) while the generating functional is independent 
of the gauge fixing. The gauge fixing also does not enter the construction of gauge invariant functions 
(such as Wilson loops). In our case, however, the gauge fixing condition is actually needed in order to 
formulate the physical time evolution and the preferred choice of gauge invariant functions on phase 
space. 

The resolution is as follows: The difference between Yang - Mills theory and generally covariant systems 
such as General Relativity that we are interested in here is indeed that the canonical Hamiltonian is 
in fact the generator of gauge transformations (spacetime diffeomorphisms) rather than physical time 
evolution. It is even constrained to vanish. In contrast, in Yang - Mills theory there is a preferred 
and gauge invariant Hamiltonian which is not constrained to vanish. In order to equip the theory with 
a notion of time we have used the relational framework discovered in [33] which consists in choosing 
fields as clocks and rods with respect to which other fields evolve. Mathematically this is equivalent 
to a choice of gauge fixing. Hence, in our case the gauge fixing plays a dual role: First, in order to 
render the generating functional less singular and secondly in order to define physical time evolution. 

2. The appearance of the 5 distributions and functional (Fadeev - Popov) determinants in (|4.17p indicates 
that we are not dealing with an ordinary Hamiltonian system but rather with a constrained system. 
One can in fact get rid of the gauge fixing condition involved if one pays a price. The price is that if 
one considers instead of q its gauge invariant extension q off the surface G = |4H I45j. then, since we 
consider the quotient Z[j]/Z[0] which leads to connected n - point functions, by the usual Fadeev - 
Popov identity that exploits gauge invariance we may replace [23J (|4.17p by 

Z\j,qi,q f ) = I [Dq Dp DX Dp] pe i s ^ x ^ e l ^ (4.22) 

Jv((ti,qiUt f ,q f )) 

However, (|4.22p is not very useful unless q(q,p) is easy to calculate which is typically not the case. 
Hence, we will refrain from doing so. Nevertheless, no matter whether one deals with (|4.17p or (|4.22p . 
the correlation functions depend on the gauge fixing G or in other words on the choice of the clocks 
[45(. [4"6] with respect to which one defines a physical reference system. 

3. The correspondence between between (|4.19p and (|4.20p also allows to reconstruct the physical inner 
product from the n - point functions: given arbitrary states tp, if)' £ TC we find a, a' E 21 such that 
||aO — V'll) ||a'f2 — are arbitrarily small. Now pick any £j < to < tf then 

< afl, a'fl >=< n, a f a'n > (4.23) 

By assumption, the operator a) a' can be written as a finite linear combination of monomials of ho- 
mogeneous degree in the components of the operator q which we write, suppressing indices for the 
components as q n . Then 

<n,q n fl>= lim <n,q(t n )..q(t 1 )n> (4.24) 

which can be expressed via (|4,19p . The existence of this coincidence limit of n - point functions is 
often problematic in background dependent Wightman QFT [22] but their existence is actually the 
starting point of canonical quantisation of background independent non - Wightman QFT as one can 
see from the identity (|4.24p . 
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4.3 The Generating Functional of Tetrad n — point functions 

We now want to apply the general framework of the previous section to General Relativity in the Hoist 
formulation. Classically its is clear that without fermions all the geometry is encoded in the co-tetrad fields 
e 1 because then the spacetime connection is just the spin connection defined by the co-tetrad (on shell). If 
fermions are coupled, the same is still true in the second order formulation so that there is no torsion. But 
even in the first order formulation with torsion one can atttribute the torsion to the fermionic degrees of 
freedom. Hence we want to consider as a complete list of configuration fields the co - tetrad. 

We will now make two assumptions about the choice of gauge fixing and the matter content of our 
system. 

I. The local measure factors p, I depends on the co - tetrad only analytically. This is actually true for 
the Hoist action [12], see also |14| . 

II. The gauge fixing condition G is independent of the co - tetrad and the Fadeev - Popov determinant 
det({-F, G}) depends only analytically on the co - tetrad. With respect to the first class Hamiltonian 
and spatial diffeomorphism constraint this can always be achieved by choosing suitable matter as a 
reference system, see e.g. |47j . However, in addition there is the Gauss - law first class constraint. 
Here it is customary to impose the time gauge gauge condition |llj which asks that certain components 
of the tetrad vanish. This will also enable one to make the connection with canonical LQG where one 
works in the time gauge in order to arrive at an SU(2) rather than G connection. Fortunately, in this 
case it is possible to explicitly construct a complete set of G - invariant functions of the tetrad, namely 
the four metricj^] g^ u = e^e^ruj and if we only consider correlators of those then we can get rid of the 
time gauge condition as indicated in the previous section (Fadeev - Popov identity). In section |4"31 
we will come back to this issue, however, in trying to make the connection of the SFM obtained with 
canonical LQG for which the time gauge is unavoidable. We will then sketch how to possibly relax 
the assumptions made under [II.]. 

Under the assumptions made we consider the generating functional xtii J) := <^]/^[0,0] where 
Z[j,J] := J [D(f> DA] 5[G[A,4>]] e^M ^(J^4>) x 

J [De] p[e,A,(t>] I[e,A,<f>] \ det[{F,G}]\[e,A,</>] e i( s ^ A ]+S m [e,AM e if M »(iAe) ( 4 . 25 ) 

Here (f> denotes the matter configuration variable. We have split the total action into the geometry (Hoist) 
part S g and a matter part S m which typically depends non trivially but analytically on e. Also the total 
current was split into pieces J, j respectively taking values in the bundles dual to those of <fi, e respectively. 

A confusing and peculiar feature of first order actions such as the Hoist or Palatini action is that from a 
Lagrangian point of view both fields e, A must be considered as configuration variables. In performing the 
Legendre transform [TJ] one discovers that there are primary constraints which relate certain combinations 
of e to the momenta conjugate to A. One can solve these constraints and then (A, e) appear as momentum 
and configuration coordinates of this partly reduced phase space. This is the reason why consider only 
correlations with respect to e. 

The idea is now as usual in path integral theory: We set 

a[e,A,cf>} :=p[e,A,<t>] I[e,A,</>] \det[{F,G}]\[e,A,<f>] e **»M.d (4.26) 

and write (|4,25p as 

Z[j, J] ■= J [D<f> DA] S[G[A, (j)}} e^M Tr (^<« <r[jjp,A, 4>] { f [De] e i s ol e A] e if M Tr(jAe)} (42?) 

11 In the presence of fermions there are additinal gauge invariant functions also involving the fermions. 
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Of course e lSm must be power expanded in a perturbation series in order to carry out the functional 
derivations with resepect to j. Indeed, if we consider just the functional integration with respect to e and 
think of A, <fi as external fields then the piece S g being quadratic in e is like a free part while S m being only 
analytic in e is like an interaction part of the action as far as the co-tetrad is concerned. Of course, in the 
computation of the physical tetrad n - point functions all the functional derivatives involved in (|4.27p are 
eventually evaluated at j = 0. 

It follows that the object of ultimate interest is the Gaussian integral 

z[j;A] := J [De]e ^ hl (4.28) 

which is computable exactly. Of course, it is not a standard Gaussian, first because the exponent is purely 
imaginary. Secondly because the "metric" Gjj(A) is indefinite so that #[j;A| would be ill defined if the 
exponent was reao. In the appendix we remind the reader how to integrate such non - standard Gaussians. 
In order to carry out this integral we must make the technical assumption that configurations A for which 
G is singular have measure zero with respect to DA. We will come back to this assumption later. 

It is at this point where we must regularise the path integral in order to perform the Gaussian integra- 
tion^! and we write the discretised version on a cubulation of M as motivated in section 14.11 that is, we 
replace (I4.28j) by the discretised version 

V,I,fJ, 

The results of appendices [A] and [B] now reveal that 

."LindfGO)) e 2 

4 : = O ' J e ~ l " (4-30) 

v vl det ( G w)l 

where we dropped a factor y^ 16 ^ for a cubulation with N vertices because it is cancelled by the same factor 
coming from the denominator in xiji J)i see (|4.25p . 



4.4 Wick Structure, Graviton Propagator and SFM Vertex Structure 
4.4.1 Wick structure 

Formula (|4.30p explicitly displays the main lesson of our investigation: The full j dependence of the gen- 
erating functional written as (|4.27p rests in (|4.30p . We are interested in the n-th functional derivatives of 
(|4.30p at j = 0. Now similar as in free field theories, the corresponding n - point functions vanish for n odd. 
However, in contrast to free field theories, for n even, the n— point functions cannot be written in terms of 
polynomials of the 2-point function. The reason is that the "covariance" G -1 [A] of the Gaussian is not a 
background structure but rather depends on the quantum field A itself that one has to integrate over. This 
renders the co - tetrad theory to be a non - quasi - free, that is, interacting theory. Nevertheless it is true 
that all Wick identities that have been derived for free field theories still hold also for the n— point tetrad 
functions albeit in the sense of expectation values or means with respect to A. 

12 As usual this prevents a "Euclidian" version of GR. Here Euclidian stands for Euclidian field theory with an analytic 
continuation to the imaginary axis of the real time variable involved (Wick rotation) which leads to a real exponent. This has 
nothing to do with Lorentzian or Euclidian signature GR. In fact, most metrics do not have an analytic section so that Wick 
rotation is ill defined and thus the connection between the real and the Euclidian theory is veiled. 

13 Actually we can formally solve the Gaussian integral w/o specifying the triangulation, i.e. we can compute it in the 
continuum. However, then one must regularise the resulting determinant which amounts to the same problem. 
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4.4.2 Graviton Propagator 

To illustrate this, consider a fictive theory in which a(e,A,<f>), G(A,<f>) are both independent of both A, e. 
This is not a very physical assumption but it serves to make some observations of general validity in a 
simplified context. This means that we can drop the <fi dependence because the generating functional 
factorises. Thus in our fictive theory we are looking at the generating functional x\j] = z \j]/ z [®] where 

r r p fiEti(G(i))) u 2 

z[j) = J [DA] z\j;A] = J J] WW")) til J |det(GW)| ] e~^^ [G ~ « (4.31) 

where is th^3 Haar measure on G. Now let 

< e£M .. e^K) >:= [ . w , m ^ ] j=0 (4.32) 

It is immediately clear that 

<<&(vi)e£(«2)>=0 (4.33) 

unless fi = t>2- This is reassuring because as we said above, physically it makes only sense to consider 
correlators of G— invariant objects such as the metric. The simplest n— point function of interest is therefore 
the 4-point function 

< 9» a v 2 (v 2 ) >=< e^(vi) ej(«i)e^(«2) e^(u 2 ) > VhJa VhJ 2 ( 4 -34) 

If we are interested in something like a graviton propagator we are interested in v\ ^ v% and obtain 

< 9^vM) W > 2 ) >= [|^] 4 < [G( Vl )-%\i\ [G(v 2 )-%lil >' (4.35) 

where for F = F[A] 

< > : = v g 4.36 

Notice that G(u) _1 does not share the symmetries of G(v), so [G~ 1 (u)]}//' ) does not vanish automatically. 
We see that we are basically interested in correlators of the inverse matrix G(w) _1 with respect to the joint 
Haar measure. Whether these have the correct behaviour in a situation where, instead of vacuum boundary 
states, one chooses coherent states peaked on a classical background metric as suggested in [381 EI] is 
currently under investigation. 

4.4.3 SFM Vertex Structure 

Finally, in order to translate (I4.36|) into spin foam language, we should perform harmonic analysis on G 
and write the integrand of the Haar measure in terms of irreducible representations of G. In particular, the 
vertex structure of a SFM is encoded in z[0] so that we are interested in harmonic analysis of the function 

fmd(G(v)) 

F(v) := ; (4.37) 
y 7 det(G(v))\ 

To derive its a graph theoretical structure it is enough to find out which F(v) depend how on a given 
holonomy A(l). Recall that F(y) is a function cylindrical over the graph 7(u) = U M<!/ 5/^(v) which is the 



14 In case of non - compact G the Haar measure is unique up to a normalisation constant which drops out in x(j). To choose 
the Haar measure instead of the Lebesgue measure makes sense in the continuum limit of infinitely "short" edges as usual. 
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union of its respective plaquette loops . Consider a fixed edge I = l^(v). It is contained in 7(1/) if and only 
if it is contained in one of the plaquette loops df^ u (v') or df Ufl (v') with \x < v or v < [i respectively. In the 
first case it must coincide either with lu{ v ') or with + v). In the second case it must coincide either 
with lu(v' + v) or with lu(v') as well. Thus in either case we must have either v' = v or v' = v — 0, v 7^ [i. 

For our illustrative purposes let us consider for simplicity that G is compact, the non compact case has the 
same SFM vertex structure but the harmonic analysis is a bit more complicated. Then each function F(v) can 
be formally expanded into SO(4) (or rather the universal cover SU(2) x SU (2)) irreducible representational 
with respect to the six plaquette holonomies A{d f ^{v)) , /i < v. These representations ir are labelled by 
pairs of half integral spin quantum numbers but we will not need this for what follows. Thus F{v ) admits 
an expansion of the form 

F ( v ) = E l 'm ■ [®m<" *r>WM v )))\ ( 4 - 38 ) 

where t'r -.is a gauge invariant intertwiner for the six - tuple of irreducible representations {t^^ v }^<v which 
is independent of v, the only v dependence rests in the holonomies. It depends on the specific algebraic form 
of F(v) which derives from the Hoist action. 

Let us define ir u ^ := ir^ for < v. By writing the six plaquette holonomies in terms of four edge 
holonomies it is not difficult to see that F{v) can also be written in the form 

F ( v ) = E L M " [®aw*" ^MW)) ® tt^A(1^v + 0))} (4.39) 

which displays explicitly the 16 variables A(lu(y)), A(l^(v + 0), v ^ \i involved and consists of 24=6 x 
4 tensor product factors. In order to arrive at (|4.39p we had to rearrange the contraction indices which 
induces the change from 1! to 1 and we have made use of ir(A(l)~ l ) = ir T (A(l)) for G = SO (4). 

We may now carry out explicitly the integrals over edge holonomies in z[0] by inserting the expansion 
P~39"|) . We write symbolical!^ 



*[0] = / II dpi H (A(l,(v))) J] F(v') 



E CQ L w] / II dMMUv))) [®v>,^u4M{iM))) ® 4M(W + m 

{ir« } v v,fi 



= EO [ / d ^9) [®^ ^0?) ® ^(ff)] (4.40) 

Here in te second step we have shifted the vertex label in one of the tensor product factors in order to bring 
out the dependence on the A(lu(v)). It follows that the end result of the integration is that for each edge 
/ = Ipiy) there is a gauge invariant intertwineJ^l 



15 This expansion would be rigorous if we knew that F(v) is an L2 function which is currently under investigation. We assume 
here that in any case we may use the Peter & Weyl theorem in a distributional sense. 

16 We rearrange the tensor products as if they were scalars but this can be corrected by performing corresponding rearrange- 
ments in the contraction structure of the intertwiners. We assume this to be done without explicitly keeping track of it because 
it does not change the vertex structure. 

17 We do not bother to expand it into a recoupling scheme and thus to label the intertwiner itself by three irreducible 
representations. 
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which intertwines six representations rather than four as in (constrained) BF theory on simplicial trian- 
gulations. The origin of this discrepancy is of course that we are using cubulations rather than simplicial 
triangulations. These six representations involved for edge l^v) correspond precisely to the six plaquette 
loops dfpviv), df^ u (v — v), v 7^ yu of which l^(v) is a segment. Therefore, if we associate to each face 
/ = ffj.u( v ) an irreducible representation 7rj = it and denote by {tt} the collection of all the 7Ty then the 
basic building block (|4.4ip can be written in the more compact form 

MM] = / dn H (g) ®icdf 7r/(ff) (4.42) 
Jg 

Likewise, if we denote ^[{-7r}] := ^{n^}i then 

*[«]=£ II ®l PllM] (4-43) 

{tt} V 

which of course hides the precise tensor product and contraction structure but is sufficient for our purposes. 

Formula (|4.43|) is precisely the general structure of a SFM. Moreover, the intertwiner (|4.42p is the direct 
analog of the intertwiner in BF theory which there defines the pentagon diagramme [6]. If we would try to 
draw a corresponding picture for our model then for each vertex v we would draw eight points, one for each 
each edge I incident at v. These edges are labelled by the intertwiner pi. Given two points corresponding 
to edges I, /' consider the unique face / that has 1,1' in its boundary. Draw a line between each such points 
and label it by 7rj. The result is the octagon diagramme, see figure CD Concretely, the edges adjacent 
to v are l/j,(v), l^(v — fi), p = 0,1,2,3. For (i / the face spanned by l^v), l u {v) is f^ u {v), the face 
spanned by lfj,(v), l v (v — v) is f^ u {v — v), the face spanned by l^(v — ft), l u (v) is f^viv — ft) and finally 
the face spanned by l^(v — ft), l v (v — v) is /^(u — fx — v). The corresponding label on the lines is thus 
n u,f} ' K ]fu 1 ' i ^v^, 7r ^u^~ u respectively. Thus the octagon diagramme has eight points and 6 x 4 = 24 lines 
(each line connects two points). These correspond to the 24 plaquettes that have a corner in v namely 
for each fi < v these are f^yiv), f^ v {v — fi), f^ u (v — v), f^v{v — fi — fi). In the case of G = SO '(4) each 
irreducible representation is labelled by two spin quantum numbers. The intertwiner freedom is labelled 
by three irreducible representations of of SO (4) and there is one irreducible representation corresponding 
to a face. Thus the octagon diagramme depends on 3 x 8 + 24=48 irreps of SO (4) or 96 spin quantum 
numbers. Since each intertwiner (|4.42j) factorises into two intertwiners [3] (one for the starting point and one 
for the beginning point of the edge but both depending on the same representations) we may actually collect 
those eight intertwiners associated to the same vertex. The collection of those eight factors is actually the 
analytic expression corresponding to the octagon diagramme which therefore maybe called the 96 j - symbol. 
The decisive difference between (constrained) BF theory and our model is however that in (constrained) 
BF theory the analog of the function F(v) is a product of 5 distributions, one for each face holonomy. 
The simplicity constraints just impose restrictions on the representations and intertwiners, but this cannot 
change the fact that there is factorisation in the face dependence. In our model, the face dependence does 
not factorise, hence, in this sense it is less local or more interacting. 

4.5 Relation between covariant and canonical connection 

Another striking feature of our model is the following: Constrained BF theory, that is, Plebanski theory, 
should be a candidate for quantum gravity. Our Hoist model should be equivalent to that theory at least 
semiclassically because morally speaking, the only difference between them lies in the technical implemen- 
tation of the simplicity constraints, modulo the caveats mentioned in section El Now one of the most 
important property of the implementation of the simplicity constraints in usual SFM is that the irreducible 
Spin(4) representations that one sums over are the simple oneJ^l. In our model we do not see any sign of 

18 If we label an irrep of Spin(4) by a pair j_) then a simple irrep. is one for which j+ = j- [6]. There is a similar 
restriction if one works with arbitrary Immirzi parameter [27| . 
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*3 




Figure 1: The octagon diagramme associated to vertex v. The eight corners correspond to the eight edges 
I = lZ{v) = ln(v + Z^-p), a = ± adjacent to v. The line between corners labelled by lZ{v), 1° (v) for \i ^ v 
corresponds to the face / = f^'(v) = f^(v + ^^fi + ^-^-u). We should colour corners by intertwiners pi 
and lines by representations 7rj but refrain from doing so in order not to clutter the diagramme. Altogether 
48 irreducible representations of Spin{A) (or 96 of SU{2)) are involved. 



that. This is an important issue because the restriction to simple representations means that the underlying 
gauge theory is roughly SU(2) rather than Spin(4) which looks correct if the SFM is to arise from canon- 
ical LQG which indeed is a SU(2) gauge theory. Thus, in usual SFM the simplicity constraints seem to 
already imply the gauge fixing of the "boost" part of the Spin (4) Gauss constraint that is necessary at the 
classical level in order to pass from the Hoist connection to the Ashtekar - Barbero - Immirzi connection 
Strictly speaking, that has not been established yet as pointed out in [38] where it is shown that the 
connection used in SFM is actually the spin connection and not the Hoist connection. But apart from that, 
in the considerations of the previous section we do not see any simplicity restrictions on the type of group 
representations. 

However, notice that what we did in the previous section was incomplete because in order to properly 
define the n - point functions we must gauge fix the generating functional with respect to the G Gauss 
constraints. Formally this is not necessary if we only consider correlators of G invariant functions such as 
the metric because the infinite gauge group volume formally cancels out in the fraction z[j]/z[0]. However, 
details matter: 

The formal arguments cannot be substantiated by hard proofs in this case. Specifically, if we consider 
G = SO(l, 3), there is no measure known for gauge theories for non compact groups (see [39] for the occurring 
complications) and thus we are forced to gauge fix at least the boost part of the Gauss constraint. This is 
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the same reason for which one uses the time gauge in the canonical theory. We expect that implementing 
the time gauge fixing [11] in a way similar to the implementation of the simplicity constraints in usual BF 
theory will effectively reduce the gauge group to SU(2). 

Details will appear elsewhere, but roughly speaking the idea is the following: 
The time gauge is a set of constraints C[e] on the co - tetrad e. By the usual manipulations we can pull the 
corresponding 5 distribution out of the cotetrad fuctional integral and formally obtain 



where XHoist is the generating functional of the previous section and xabi stands for the Ashtekar - Barbero 
- Immirzi path integral. 

Whether this really works in a rigorous fashion remains to be seen. However, we find it puzzling that 
the simplicity constraints in usual SFM, which classically have nothing to do with the time gauge, should 
automatically yield the correct boundary Hilbert space. It seems intuitively clear that the time gauge must 
be imposed in the quantum theory in addition to the simplicity constraints, just like in the classical theory, 
as we suggest. Without imposing it, we do not see any sign of a restriction from G to SU(2) in our model 
where we solve the simplicity constraints differently. We expect this to be related to the work [50]. This 
observation indicates that usual SFM and our formulation are rather different from each other. 

5 Conclusions and Future Work 

In this paper we have proposed a different strategy to construct spin foam models for GR. Rather than the 
Plebanski action we take the Hoist action as our starting point. This means that the simplicity constraints 
of the Plebanski formulation have been correctly taken care of. The price to pay is that the connection to 
BF theory is lost. The motivation behind our strategy is that BF theory is a TQFT and therefore quite 
different from GR which has an infinite number of physical degrees of freedom. Hence the usual triangulation 
independent methods developed for TQFT and employed in current SFM are possibly less powerful in the 
context of GR. In particular, the fact that it is difficult to deal with the simplicity constraints in current 
SFM might be a sign of that. Another problem with the Plebanski formulation that we have not mentioned 
yet is that it is difficult to couple matter because matter directly couples to the cotetrad rather than the 
B field. In principle one can express e via B modulo simplicity constraints but the corresponding formulas 
are even more involved than those for e. Notice that one must couple matter to BF theory in order to get a 
realistic model. For 3D gravity coupling matter is straightforward [51] because there B field and e coincide 
while in 4D this has not been done yet except for non standard model fermions which just couple to the 
connection [52] or membranes coupled to pure BF theory [53] . 

The method we proposed in this paper might be a called a brute force and textbook strategy. Dropping 
the insistance on triangulation independence right from the beginning we proposed a Wilson action like 
naive discretisation of the Hoist action. We carefully studied the connection of the Hoist path integral with 
the canonical LQG correlation or n - point functions and used relational techniques to make the connection 
with the physical Hilbert space and observables. In principle, none of these ingredients are new, they have 
been successfully employed in other contexts. Of course, the appealing elegance of (constrained) BF theory 
has disappeared in our formulation, the integrals to be computed are rather challenging (but not impossible) 
and the gauge fixing conditions for spatial diffeomorphism and Hamiltonian constraint as well as a local 
measure factor without which the connection to observation is lost complicate the formalism. Yet, we feel 
that the resulting structure, while far from being worked out in detail, has some interesting features such 
as the Wick like structure of the physical tetrad correlators and a less local structure. In particular, we 
have shown that imposing the time gauge also in the quantum theory comes out as a necessary and natural 
condition in our model in order to make contact with the LQG Hilbert space. 

As already mentioned in the introduction, this paper is exploratory in nature. It focuses more on ideas 
rather than analysis and there are many open issues that need to be settled before the present model can be 




(4.44) 
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taken seriously. Apart from the topological issues mentioned in section 14.11 the convergence and measure 
theoretic issues discussed in section I4T31 and finally the issues with the imposition of the time gauge outlined 
in section [431 there are further points that need to be addressed. One of the most serious ones is the contin- 
uum limit: The fact that we are working with cubulations suggests a naive but natural notion of continuum 
limit which consists in studying the behaviour of the correlation functions under barycentric refinement of 
the hypercubes at fixed IR regulator (boundary surface). Of course, in the spirit of the AQG framework 
|20| one could also say that the continuum limit has been taken already provided that one works with 
infinite cubulations. This requires then, in a seperate step, to remove the IR regulator. A more practical 
but still important problem is the following: Even at finite UV and IR regulator, it is already hard enough 
to compute the determinant of the covariance matrix of the co-tetrad Gaussian and to determine its index 
(which may vanish automatically, see appendix \B§ . But since these covariances are highly correlated, the 
practical computation of the n - point functions at least in the macroscopic regime will be possible only if 
the corresponding non trivial measure has some kind of cluster property [S3]. We hope to come back to all 
of these issues and many more in several future publications. 
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A Non — standard Gaussian integrals 

Let G be a real valued, symmetric and non singular matrix on the real vector space V = M. n and let j £ V. 
We are interested in the non - standard Gaussian integral 



/ := / d n x e% x Gx+j x] (A.l) 
Jv 

Using the translation x = y — G~ 1 j we can simplify this to 



There exists an element S E GL(n, K) such that G = S T D S where D is a regular diagonal n x n matrix 
which is possibly indefinite. Denote by d\, .., d n 6 M — {0} the entries of D. Then the change of coordinates 
y = Sz reveals 

' ' d n zet z * Uz (A.3) 



|det(5)| 

Now consider the basic integral 

I d := [ dze ldz2/2 (A.4) 



V 



for d £ M — {0}. For d = ik, k > we know the value of (|A.4p . however, that formula involves a square 
root and thus analytic continuation of Id in d is ambiguous. Hence we must determine the value of (|A.4p 
by independent means. 

The integrand in (|A.4p is entire analytic in z without poles. For d > or d < respectively the integral 
over the arc < arg(z) < ir/2 or > arg(z) > —ir/2 respectively vanishes at infinite radius. Hence, 
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using a Cauchy integral argument we may rotate the integral from z E K to z G e ,sgn ( rf ) 7r 
z = e^gn^W^ f G R we get 



so that with 



(A.5) 



Given a symmetric matrix G with signature q (i.e. p positive, g negative and n 
we define its index ind(G) := p — q. Then, combining (|A.3P and (|A.5P we obtain 



2vr e 



-i|i T G- 1 j «nd(G)7r/4 



p — g zero eigenvalues) 



(A.6) 



|det(G)| 

B On the index of special matrices 

While determinants maybe tedious to calculate, its is always analytically possible. However, the index 
of a matrix is harder to obtain. While there exist algorithms to obtain it just from its characteristic 
polynomial (rather than from its spectrum which would be impossible to determine analytically for general 
large matrices) for concrete matrices such as Descartes sign rule |55| , for general matrices of a given restricted 
structure there are no such algorithms available except for in a few cases. 

Our situation is the following: Consider the 16 x 16 matrix G with entries Gr^j\r v j\ := G^j. Since 
G v j\ it is symmetric. Let us also write := el. We consider the lexicographic 
ordering of the compound index (jjj) as (00), .., (03), (10), .., (13), (20), .., (23), (30), .., (33). Consider the 
antisymmetric 4x4 matrix G^ u with < fi < v < 3 given by {G tiV )ij := G^j. Then the 16 x 16 matrix G 
has the following block structure 



G 



This leads us to the following conjecture. 
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Conjecture B.l. 

Let A, B, C, D, E, F be antisymmetric, real valued 4 x 4 matrices and let G be the symmetric 16 x 16 matrix 
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D 
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-B 


-D 
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-C 


-E 


—F 


o ) 



(B.2) 



Then ind(G) = 0. 

It turns out to be extremely hard to prove this conjecture although it is rather plausible. For example, 
it is easy to show that the conjecture is correct when the matrices A, B, C, D, E, F are 2x2 antisymmetric 
matrices. It is also true when the matrices A, B, C, D, E, F are linearly dependent. We delay the proof (or 
disproof) of this conjecture to future publications. 

If the conjecture was true and G is non singular then we would know that G has eight positive and 
eight negative eigenvalues. Hence we would know that det(G) > 0. In order to compute det(G) we make 
use of the following basic factorisation property for an arbitrary block matrix with blocks A, B, C, D with 
invertible A 

A B \ _ f A W 1 A- X B 
C D ) V C 1 ) [ D — CA- l B 



(B.3) 
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It follows that 

det(G) = det(A) det(D - CA^B) (B.4) 

In our situation, by means of (jB.4[) we can iteratively downsize the size of the matrix of which we have to 
compute the determinant from rank 16 to 8 and then to 4. At rank 4 we may use Cayley's theorem [55] 
in order to express det(G) directly in terns of polynomials of the traces of products of the G^ v and thus in 
terms of traces of products of the plaquette loops. 
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